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ABSTRACT 

With abundances and 2-cr upper limits of fluorine (F) in seven metal-poor field 
giants, nucleosynthesis of stellar F at low metallicity is discussed. The measurements 
are derived from the HF(l-O) R9 line at 23358 A using near-infrared ^-band high- 
resolution spectra obtained with CRIRES at the Very Large Telescope. The sample 
reaches lower metallicities than previous studies on F of field giants, ranging from 
[Fe/H] = -1.56 down to -2.13. Effects of three-dimensional model atmospheres on 
the derived F and O abundances are quantitatively estimated and shown to be insignif- 
icant for the program stars. The observed F yield in the form of [F/0] is compared 
with two sets of Galactic chemical evolution models, which quantitatively demon- 
strate the contribution of Type II supernova (SN II) v-process and asymptotic giant 
branch/Wolf CRayet stars. It is found that at this low-metallicity region, models can- 
not well predict the observed distribution of [F/0], while the observations are better 
fit by models considering an SN II v-process with a neutrino energy of £'v = 3 x 10^^ 
erg. Our sample contains HD 110281, a retrograde orbiting low-or halo star, showing 
a similar F evolution as globular clusters. This supports the theory that such halo stars 
are possibly accreted from dwarf galaxy progenitors of globular clusters in the halo. 

Subject headings: Galaxy: halo — stars: abundances — stars: Population n 
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1. Introduction 



Although fluorine is a very interesting element, its chemical evolution is not yet well 
understood. F has only one stable, yet fragile isotope, ^^F, which is easily destroyed in the 
stellar interiors by the most abundant elements, H and He, by means of the reactions '^F(p,Q')'^0 
and ^^F(Q',p)^^Ne. Theoretically, in order to explain the production of F, three sites have been 
proposed to explain the mechanism that enables it to escape from the hot stellar interior after it 
has been created, including Type II supemovae (SN II), asymptotic giant branch (AGB) stars, and 



spallation ( 


Wooslev & Haxton 


1988; 


Wooslev & Weaver 


1995 


Heaer et al. 


2005 


Kobavashi et al. 


2011a 


). For AGB stars. 


Forestini et al. 


(1992 


) proposed that ^^F is synthesized and then dredged 



up to the surface during the He-burning thermal pulses. In WR stars, '^F is probably produced 
during the He-burning phase and th en ejected into space by strong stellar winds before it is 



destroyed (jMeynet & Amould 



2000) 



There are very few atomic or molecular lines of F that are suitable for stellar abundance 
analysis. Current studies on F have focused o n the measurement of the single HF molecular line 



at 23358 A in cool stars. Jorissen et al. 



(| 19921) put up a pioneering work in the field by measuring 
the first F abundance of a number of giants with near-solar metallicities. Through the past decade, 
more studies on F nucleosynthesis have been carried out, aimi ng at various types of objects. 



including solar- and mo derate-metallicity global cluste r giants (ICunha et a 



2005 



Smith et al. 



2008). AGB stars (Abia et al. 



2005h. giants in the G alactic bulge (|Cunha et al 



(Recio-Blanco et al. 



2009 



2008 



2003 



Wemer et al. 



Uttenthaler et al. 



2010|) . and main-sequence dwarfs in the solar neighborhood 



2012h . Besides, detailed exploration of the chemical evoluti on of F 



las been 



carried out targeti ng chemically 



C-enhanced stars (Schuler et al 



pecu 



2007 



and even an extragalactic C-star (|Abia et al 



i ar objects, such as co ol extreme He stars (iPandev 



Lucatello et al. 



2011 ). 



201 ih . Ba stars (lAlves-Brito et al. 



2006), 



20111), 
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With new nucleosynthesis yields and updated observations, chemical evolution models 
can now predict more details of the nucleosynthesis of F, but only the contribution of AGB 



stars t o the F production has been observationally confirmed dWemer et al 



2005 



2010, 



Abia et al. 



Renda et al. 



2011|) . The chemical evolution model for MilkyWay-like disk galaxies of 
(|2004|) suggests a scenario where AGB stars only began to contribute to the production of F 
after ~ 1 Gyr. This model states that AGB stars are at present the major contributors for the 
production of F; while in the early Galaxy, only SN II contribute to the F abundance, and provide 
a solution to the discrepancy between models and observed F abundance in Galactic giant in the 
intermediate-me t allicity range. A more recent chemical evolution model targeting fluorine by 



Kobavashi et al 



(1201 lal ) adopts the v-pro cess yields and yields of A GB stars with 1-7 Mq, based 



on the infall and star formation history in 



Kobavashi et al. 



(1201 Ibh . The model predicts a [F/0] 



plateau at low metallicity depending on the neutrino luminosity. The included v-process enhances 
the predicted F abundances at higher metallicities, and helps to fit the high F abundances observed 
in field and bulge giants. Both models are in good agreement with observations in regions with 
solar or moderately low metallicities, indicating the importance for c onstraining the theoretica l 



Alyes-Brito et al. 



( 20121) 



yields at low metallicity for further scrutiny and future observations, 
filled the gap at the low-metallicity area with M22 data, but the complex chemical enrichment 
history of this globular cluster leads to a large dispersion (0.6 dex) of the F abundances, thus does 
not follow the trend of other globular clusters or being explained by a general model. Hence, 
the metal-poor field giants in our observed sample, which have metallicities of [Fe/H] < -1.5, 
become very interesting in allowing us to map the pristine F abundance in earlier phases of the 
Galactic halo evolution. 

Observations and data reduction, including the determination of stellar parameters and 
abundance analysis, are addressed in Sections |2]through |4] The three-dimensional (3D) eff"ects on 
the abundances are described in Section |5] In Section |6l our results are discussed and interpreted, 
and a brief summary is given in Section |7] 
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2. CRIRES Observations 

Seven metal-poor field gia nts were observed with CRIRES P (Cryogenic high-resolution 



Infrared Echelle Spectrograph; 



Kaeufl et al. 



20041) at Very Large Telescope UTl in a number 



of service runs from 2008 May through 2009 March. The slit width was set to 0".2, yielding 
the maximum resolution of 100,000 (~ 3 km s~'). The four CCD chips cover the wavelength 
ranges from 22993 to 23118 A, 23152 to 23273 A, 23304 to 23420 A, and 23450 to 23560 A. 
Chip 3 includes the HF(l-O) R9 line which, in cool stars, has been verified to be most reliable 
for F abundance determination. The integration time was 30 s or 45 s for each nodding position 
per setting, and for two of the program stars, two exposures were obtained. A hot standard star 
with similar air mass was observed immediately either beforehand or afterward. The raw frames 
of science and standard objects were reduced with the CRIRES pipeline (v2.0.0), and after that, 
the science spectra were divided by corresponding standard star spectra to correct for the telluric 
lines and the illumination pattem, using the IRAF task Telluric. The typical signal-to-noise (S/N) 
ratio of the sample spectra around 23358 A is higher than 200. Spectra of the HE line area for all 
exposures of our program stars are shown in Figure [1] 



3. Atmospheric Parameters 

The eff'ective temperature T^g and the surface gravity log g can be determined either 
spectroscopically or photometrically. The microturbulent velocity ^ is normally derived by 
removing any trend of abundances with equivalent widths of Fe lines. However, the fact that 
there is no Fe feature in the sample spectra within the coverage of CRIRES observation makes it 
difficult to independently determine the atmospheric parameters with the IR spectra alone. Thus, 
the stellar parameters for the sample stars were adopted from previously published data. Out of 



^ http : w w w. eso . org/instruments/crires 
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a number of measurements, we relied on the work of IShetronel (|l 996|) . which includes all seven 
objects in our sample and thus provides a uniform set of atmospheric parameters. Moreover, it 
contains measurements of equivalent widths. 



The model atmospheres adopted by 



Shetrone 



(|l996h are taken from the grid of 



Bell et al 



(119761). which i s now outdated. Therefore, we have adopted the atmospheric parameters of 



Shetronel (|l996|) as the initial values, and further adjusted ther n with local th ermal equilibrium 



1994 . The adjustment 



(LTE) one-dimensional (ID) models computed with ATLAS9 (iKuruczl 
followed the traditional spectroscopic method. Tes was adjusted to make Fe I lines with low 
excitation potential yield the same abundances as those with high excitation potential; the 
microturbulent velocity was adjusted until strong Fe I lines yield same abundances as weak 
Fe I lines; and log g was adjusted so that abundances from Fe I and Fe II lines were forced to 
agree with each o ther. The equiva lent width (EW) of the Fe lines were directly adopt ed from th e 



measurements of 



Shetronel (|l996|) . The adjusted parameters do not differ much from 



Shetrone 



(Il996h 's values, with a difference being 10+45 K for T^s, -0.03+0.08 for log g and -0.02+0.03 



for [Fe/H]. 



Shetronel (Il996[) . abundances of 



Similarly, with the EW measurement and the linelist of 
O, Mg, Al, and Eu (which cannot be determined from CRIRES spectra since no corresponding 
line features are present in the ^-band spectra) were also re-calculated, using the ATLAS 9 
atmospheric models. The re-derived abundances show overall agreement with Shetrone's values, 
e.g., with a difference of 0.10+0.08 for O, 0.10+0.06 for Mg, -0.09+0.07 for Al, and 0.06+0.06 
for Eu, respectively. The adopted atmospheric parameters and abundances, based on the optical 
data, are listed in Table [Tj 
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4. Abundance Analysis 



with ATLAS9 (Kurucz 



Using the spectrum synt hesis code MOOG (ISnedenlll973|) and model atmospheres calculated 



19941) . LTE Na, C, and F abundances (and upper limits) were obtained 
by comparison with CRIRES spectra. The mixing-length parameter was set to aMLT = 1-25, and 
no convective overshooting was adopted. The wavelength range of the CRIRES spectra covers 
a few Phillips system C2 lines; however, since our sample stars are of rather low-metallicities 
and without carbon enhancement (except for HD 135148), no C2 line was detected. Therefore, 
we relied on the O abundances obtained as previously described, and derived the C abundances 
by fitting several fi rst overtone vibrational-rotational lines (2-0 and 3-1 vibration series). 



Although iShetronel (|1996h also provides measurements of Na abundances, which was determined 
as an average of the abundances derived from EWs for the 6154 A and 6160 A Nal lines and from 
synthetic spectrum fits to the 5682 A and 5688 A Nal lines, there are two objects without Na 
EW measurements in the literature. However, there are two Na lines located at 23379 A, which 
are strong and clean enough to measure the Na abundance with the IR spectra. The left plot of 
Figure [2] compares Na abundances determined from CRIRES spectra and re-derived from the 



EWs measured by 



Shetronel (|l996h for the five stars with EW measurements, an d the righ t 



plot 



compares the CRIRES Na abundances with the average Na abundances given by Shetrone 
The measurements agree within their uncertainties, and the CRIRES values are adopted for the 
following discussion. 

The F abundances of our sample stars were determined by fitting the unblended HE (1-0) R9 
line at 23358 A. Out of the seven sample objects, clear HE R9 features were detected in two of 
them, HD 1 10281 and HD 135148, as shown in Figure[3l Note that HD 135148 i s a spectrosco pic 



binary with a companion star with the approximate mass of a white dwarf (e.g., Carney et al. 



20031), it is thus not likely that the detected HE feature is reflecting the closeby CO lines from 
the companion stellar spectra. To confirm this, we have also checked the auto-correlation of the 
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spectra, and found no significant second peak associated with a companion star. For the remaining 
five spectra, 1-cr upper limits of F abundances were obtained, adopting the statistical error of EW 



based on the classical formula of 



Cavrel 



( 1988 ) 



{t^Wy^ ^ l.6{w^xf'^e (1) 

where w is the FWHM of the line; Ax refers to the sampling step of the spectra, and e to the 
reciprocal spectral S/N in the case of normalized spectra. Considering the fact that the HF line 
is weak in the CRIRES spectra, we can estimate their expected FWHM from weak CO lines 
close to the HF line. The S/N of the spectra were calculated from several continuum regions 
on the same CCD as the HF feature. Note that except for the spectroscopic binary HD 135148, 
two other objects, HD 003008 and HD 1 10281 clearl y show broad e r line s than other program 



Carney et al. 



( 20031) have concluded that 



stars. Spectroscopic and photometric observations by 
the line-broadening of these two stars is caused by fast rotation rather than stellar companions, 
and also proposed the idea that these two rapid rotators show enhanced mass loss which is 
aided by the rotation, or by the addition of both orbital energy and angular momentum from an 
accreted planet to the stellar envelope. We have also investigated the detectability of HF lines with 
various EWs/upper limits, e.g., the difference in A(F) between the two exposures of HD 1 10281 
corresponds to a 1-2 mA difference in EW, which is on the same level as the 1-cr upper limit. We 
have adopted 2-cr upper limits of A(F) in the following discussion, to ensure the possible detection 
considering the typical S/N and line widths of our program stars. 

For the only C-enhanced star HD 135148, the ratio of '^O/^^O was derived based on three 
i^C^'^O lines in our spectra, at 23347 A, 23357 A, and 23379 A. A ^^O/^'^O ratio of 140+50 was 
derived. This oxygen isotope ratio is higher than normal at such a metallicity, if one only considers 
the nucleosynthesis of massive stars. However, since '^O can be overproduced by AGB stars, 
and AGB stars can be the typical companion star that transfers nucleosynthesis products to the 
C-enhanced star, the high ratio of '^O/'^O can be explained by this scenario. Note that the '^C'''0 
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line at 23357 A is located right to the blue wing of the HF R9 l ine, bu t is not at all affecting the 



derived F abundance, as testified in early work like 



Abiaetal 



(120091) ■ The fitted ^^C^^ Q line is 



also s 



10 wn in Figure |2] For the synthet i c fittin g, we have adopted the N a linelist from 



(119941) . the CO linelist from 



GoorvitchI (Il994h . and the HF linelist from 



Cunha et al. 



Kurucz 



(2003). 



Validity of the line list has been checked by fitting the ^-band spectra of a Boo from the atlas of 



Hinkle et al. 



( 1995 ) 



The uncertainties of the abundance of F and other elements mainly comes from two aspects: 
the statistical error in the line equivalent widths, and the uncertainties of the atmospheric 
parameters. To estimate the former, statistical errors of the EW measurements were calculated 
adopting Equation [H and the latter were determined by individually varying stellar parameters 
to measure their separate effects on the derived abundances. Typically, e.g., for HD 110281, 
the 1-cr EW error leads to an uncertainty of 0.02 dex of the F abundance, AT^s = +100 K leads 
to an increase of +0.25 of the F abundance, Alog^ = +0.5 produces AA(F) = -0.02, and 
= +0.5 km s"' gives a negligible 6A{F) < 0.001. As for the Na abundance, the typical 
sensitivities to EW errors and atmospheric parameters are AA(A^a) = +0.06, +0.07, -0.03, and 
-0.01 respectively. The resulting uncertainties are then summed in quadrature to obtain the 
total statistical errors. For stars with more than one exposure, i.e., HD 110281 and [S84] 2643, 
abundances for the exposure with higher S/N was adopted, and the differences between the derived 
element abundances from different exposures turned out to be negligible. The derived abundances 
of C, O, F, Na, Mg, Al, and Eu are listed in Tabled 
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5. 3D effects 



5.1. Abundance Corrections for the HF R9 and Forbidden O Lines 



It is well-known that atmospheres of metal- depleted, late-type s tars are prone to exhibit 



significant deviations from radiative equilibrium ([Asplund et al. 



1999|). In particular, abundances 



Gonzalez Hernandez et al. 



deriv ed from molecular species can experience large corrections (e.g., 

lOOsi . From this perspective it is desirable to quantify the 3D effects - at least for fluorine and 
oxygen. Unfortunately, we have no 3D models available at the rather low gravities of our program 
stars. To nevertheless derive an estimate of the abundance corrections, we took four 3D models 



from the CIFIST atmosphere grid ([Ludwig et al. 



2009|) at higher gravity and inspected the trend 



of the abundance corrections with changing T^s and lo gg. When ca ' 



Caffau et al 



culati ng the abundance 



(1201 II) . In the 3D models 



corrections, we followed the methodology described in 
we assumed a solar chemical composition scaled to a metallicity of [M/H] = -2.0, with an 
enhancement of the or-elements by -1-0.4 dex. 

Table |2] summarizes the results for the HF R9 line, and the two forbidden oxygen lines at 
6300 A and 6363 A. The listed 3D abundance corrections are evaluated for weak lines (negligible 
saturation) - a good approximation considering the typical line strength. When running the 3D 
models, sufficient statistics were gathered so that residual statistical uncertainties on the abundance 
corrections are smaller than the accuracy to which they are given in the table. For completeness, 
we mention that in the ID hydrostatic comparison models we assumed a mixing-length parameter 
of ttMLT =1-0 (close to what was adopted in the ATLAS9 models), and in the related ID spectral 
synthesis calculations a microturbulence of ^ = 2.0km s However, both aspects have hardly 
any impact on the result. The corrections in Tab. |2] suggest that for the low-gravity giants of our 
sample, minor 3D effects are to be expected. Only at gravities log ^ > 2 we expect substantial 
corrections (> 0.1 dex) on the F abundance at [M/H] = -2. For the HF line, mainly horizontal 
inhomogeneities (as opposed to differences in the vertical stratification) are the cause of the 
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changes of its strength in 3D compared to ID models. 



Digressing for the moment, we note that the abundance correction s obtained for the fo rbidden 



Collet et al. 



( 2007 ). These 



oxygen lines are not fully compatible with the corrections obtained by 
authors obtained ^ -0.1 dex at T^ff = 5050 K, logg = 2.2, [M/H] = -2.0, and [0/Fe] = +0.5. 
Even considering that the atmospheric parameters are not exactly the same in their and our 
models, the correctio ns of Collet and collaborators appear noticeably la rger. This is also the case 



for other species (see 
understood. 



Ivanauskas et al. 



20101 : 



Dobrovolskas et al 



2010|) and the reason is not yet 



Nevertheless, we proceed under the assumption that 3D corrections are insignificant for the F 
and O abundances we derived for our program stars. 



5.2. Convective Line Shifts of HF and CO Lines 

We used the cores of strong CO lines as synthesized by a ID spectral synthesis code to 
correct for the imperfections of the wavelength calibration provided by the CRIRES reduction 
pipeline. One might ask whether convective line shifts alter the wavelength position of the cores 
significantly with respect to their laboratory wavelength. Again, we do not have fully adequate 
3D models available for our target stars. We nevertheless used the 3D model from the CIFIST 
grid which gets closest (Teff = 3886 K, log^ = 1.0, [M/H] = -2.0) to calculate line bisectors of 
the HF R9 line and a prototypical CO line in the HF region. We varied the strength of the lines by 
(artificially) varying the oscillator strength; in the case of the CO line this means that the variation 
was done at fixed C/0 ratio. Figures |4]and|5]illustrate the outcome. The perhaps surprising result 
is that the magnitude of the convective shifts is very modest, not exceeding 250m s The cores 
of strong CO lines in our target stars have typically a residual flux around 0.7, i.e., the HF line 
is weak. Diff"erential shifts between their cores are even significantly less, unlikely to exceed 
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1 ° 

100 m s" corresponding to « 8 mA at the wavelength of the CRIRES spectra - much smaller than 
the line width. We conclude that our wavelength correction procedure should produce reliable line 
positions for the HF and CO lines, as can be seen in Figured! 



6. Discussion 



With the measurements of F abundances and upper limits of the program stars, it is now 
possible to probe the behavior of F evolution of field giants with low metallicities. Figure |6] 
displays the kinematic properties of our program stars. Please note that the radial velocity and 
proper motion data are directly obtained from previous literatures as listed in the database of 
SIMBADy. As can be seen, one of the stars with clear HF detections, HD 110281, is identified 
as a halo star with retrograde rotation, and another program star, BD +01 2916, is classified 
as a thick disk star The sensitivities of the F abundances are checked along temperatures, and 
abundances of C, O, and Fe, as shown in Figure |7] Relevant c orrelations or dependencies have 



been addressed in a nu mber of articles, including observations (|Jorissen et al 



1992 



Cunha et al 



2003 



Abia et a 



models such as 



2009h. as well as the oretical predictions frorn Galactic chemical evolution 



Renda et al. 



(2004) and 



Kobavashi et al. 



(1201 lal) . We would expect for giant stars 



to find a correlation between the abundance of F and the abundances of O and Fe, and that the F 
abundances could be enhanced by stars exhibiting C-enhancement. It is quite likely that potential 
trends as expected are present in the comparisons, but we cannot draw firm conclusions on this 
since our sample size is small and we have obtained mostly upper limits only. Nevertheless, it is 
clear that HD 1 10281, the object with the highest F abundance, is the star with the lowest Teff in 
our sample, and showing the highest metallicity (Fe and O). The HF molecule formation is very 
sensitive to temperature, hence T^s is expected to set the level of upper limits on the F abundance. 



^http : //simbad . u- strasbg . fr/simb ad/ 
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The other star with a HF-detection, HD 135148, is showing notabl y higher C abundan ces 



compared to o t her p rogram stars, since it is a CH-star as shown by 



Shetrone et al 



Carney et al, 



( 2003 ) and 



(|1999|) . The observed trend of the F abundances with metallicity tracers suggest 
that O and Fe may be used as a reference on the detection limit of F; however, it is not clear to 
us whether there i s a low er limit of F detection determined by the stellar metallicity, because 



Alves-Brito et al. 



(|2012|) has studied several M22 giants with metallicities (and stellar parameters) 



that are comparable to our sample but most of their giants' spectra s how measurable HF features. 



Note that these spectra based on which the 



Alves-Brito et al. 



(120121) derive the F abundances do 



not present higher resolution or S/N than our sample spectra; therefore, it is suspected that the 
different detection fraction of HF features between the two samples are possibly due to the fact 
that one is in clusters while the other is in the field. 

In Figure [H [F/0] abundance ratios treated as the indicator of the F yield of the sample, are 
plotted against [0/H] and [0/Fe] and compared with model predictions and previous observations. 
Two sets of models are selected, including Renda et al. (2004, R04 hereafter) which for the first 
time quantitatively demonstrated the contribution of AGB and WR stars to the F production, and 
the latest model by Kobayashi et al. (201 la, hereafter Kl 1), who calculate the contribution of 
the v-process of SN II and hypemovae (HN) to the nucleosynthesis of fluorine. Three model 
variants are constructed in R04: MWa (thick solid line) with SN II as the only source of ^^F, MWb 
(thick dotted line) including yields from both SN II and WR stars, and MWc (thick dashed line) 
including all three sources of SN II, WR, and AGB stars. For all three cases, R04 have considered 
that neut rino spallation is largely do minant in fluorine production, and they adopted the SN II 



yields of 



Woosley & Weaverl (|1995|) . which assumes a total energy in neutrinos of Ey = 3 x 10^^ 



erg. As for Kll, five cases are included: evolution with neither AGBs nor v-process (SN Il/Ia 
-I- HN, triple-dot-dashed line), chemical evolution with AGB contributions (SN Il/Ia -i- HN -l- 
AGB, dashed line), the model for globular clusters with v-process (dash-dotted line), v-process 
of SN II and HN (dash-dotted and solid lines) with neutrino luminosity of Ey = 3 x 10^^ erg and 
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Ey = 9 X 10 erg respectively. Note, that Kl 1 models do not include yields of WR stars, because 
the contribution by WR to F may be reduced by including rotation in the stellar models. For more 
details of the models, we refer the interested readers to R04 and Kl 1. 

As can be seen in Figure [81 for all cases taking v-process into account, all models predict 
a plateau of [F/0] at low metallicity with [0/H] < -1.2, as proposed by Kll. The v-process 
remarkably enhances the level of F production compared to models without y-process and 
indicates the dominant contribution from neutrino spallation to producing F in the early stage of 
Galactic evolution. For all four model variants with the same neutrino energies of fiy = 3 x 10^^ 
erg, the predictions of the F production at the plateau of the models of R04 and Kll agree 
well with each other. For a better judgment of the significance of the differences between our 
observations and the model predictions, 1-cr error bars are also included to the upper limits of 
[F/0] of the sample. As C-enhanced stars are believed to be s-process enhanced and producers 
of F, the CH-star HD 135148 shows the highest F abundances as shown in Table [H thus is 
not included in Figure [8] which intends to compare with the general trend. Also noted, due to 
its different kinematical properties from the rest of the sample, HD 110281 (shown in filled 
circle) may not follow the Galactic halo evolution, hence should be separately considered when 
compared with models. Except for these two objects, all upper limits of the other five objects in 
our sample (as shown with lower triangular) are located around the predicted plateau of R04 and 
the SN II v-process model of Kl 1, while well below the HN model of Kl 1 with = 9 x 10^^ erg. 
The general trend of the observed [F/0] with increasing metallicity, and the apparent discrepancy 
between the observed [F/0] and predictions by models without neutrino spallation, suggest that a 
contribution of the v-process is needed to explain the early production of F in the Galactic halo. 
The program stars with intermediate F content in the sample tend to follow the evolution model 
for SN II v-process of Kl 1 and the early phase of R04 models, which predict a similar "plateau", 
but program stars with low F upper limits are distributed at relatively lower level of [F/0], and 
cannot be explained by these models. 
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Our observed F abundances are also compared with previous observations, including gian ts 



of glo bular clusters (ICunha et al. 



2003 



2012 ) , the Milky Way 



2003 



Smith et al. 



Cunha & Smith 



bulge (iCunha et al. 



2005 



Yong et al. 



2008 



Alves-Brito et al 



20081), as well as K and M field giants (|Cunha et al 



2005b . As show n in the left panel of Fig ure El for globular cluster giants 



with comparable metallicity, e.g., M22 (jAlves-Brito et al 



2012h . the detected F level shows clear 



diff'erence from our program stars, which could be explained by the different chemical enrichment 
history of field and cluster stars, and/or the wide range of chemical abundance variations of 
M22 due to its complex formation history. Previously available data of field giants covered only 
metallicities close to solar, thus they show a notably higher F abundances than our metal-deficient 
field giants. Apart from the field and globular cluster halo giants. Milky Way bulge giants 



(ICunha et al. 



2008h are also compared. They are showing remarkably higher F abundances 



with a large scatter, which is difficult to explain by yields of current models. This also implies 
a diff'erent evolution of this central component of the Milky Way. However, it is interestingly 
noticed that HD 1 10281 follows very similar trend of [F/0] as giants in globular clusters such as 
M4, NGC6712, and to Cen. Note in Table [U that this object also shows rather low Na and Mg 
abundances, as w ell as being on a retrograde orbit in Figure |6l thus it is very likely a "low-a" halo 



star as defined by 



Nissen & Schuster (2010 



201lh . As they have found based on the abundance 



and kinematic analysis, the less-bound low-or stars are probably accreted from dwarf galaxies with 
relatively slow star-formation, and some of them may be associated with the progenitor galaxy of 
globular clusters, including to Cen, which could refer to a similar or relevant chemical evolution, 
hence explain the similarity in F evolution of HD 110281 and the compared globular cluster 
giants. Future observations of abundances of F and other elements of similar objects would be 
very interesting and important to obtain a deeper understanding of the origin of the dual halo of 
the Galaxy. 

The timescale of SN II is much shorter than that of AGE or WR stars, which naturally 
explains that SN II dominate the area with lower [F/0] and higher [0/Fe], whereas AGB stars 
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start to notably contribute after [O/H] > -1.5 and become more important at higher metallicity 
and higher [F/0]. Therefore, the F abundance is a good clock to distinguish the contribution from 
low-mass AGB stars and SN II. It is true that based on our comparisons, none of the included 
models could well fit or predict the observations at lower metallicities, neither for field giants nor 
globular clusters. But it is also true that both the R04 and Kl 1 models are originally aimed at 
better fitting the observational data of globular clusters at intermediate metallicities and giants 
toward solar metallicity due to the availability of observational data. This is also the reason why 
the observed F production from AGB stars is better explained, while there is still little known 
about the early evolution of F dominated by massive stars or the different F evolution of halo field 
giants and globular clusters stars. Future observations of field stars at low metallicities are crucial 
for further constraining models for such early stages of the Galactic chemical evolution, e.g., 
the neutrino luminosity released from SN II, and the initial conditions of the chemical evolution 
models for the older halo. 



7. Summary and conclusion 

With CRIRES observations of seven halo field giants at [Fe/H] < -1.5 and [O/H] < -1.1, 
and abundances of O, Mg, Al, and Eu determined in optical spectra from the literature, we have 
investigated the nucleosynthesis of their F content, and for the first time, explored the evolution 
of F in field giants in the low-metallicity regime. Among the program stars, the HF R9 feature 
is detected in two cases, HD 1 10281 and HD 135148, and 2-cr upper limits are estimated for the 
F abundance of the other stars. Using CIFIST model atmospheres, the 3D eff'ect on the F and O 
abundance with comparable stellar parameters to the program stars are estimated, showing that 
for our sample, the corrections are negligible. The sensitivity of the F abundance/upper limits are 
also checked against T^s, metallicities, O, and C. A correlation can be found between the upper 
limits of F and T^g, which confirms the expectation that more stringent limits can be obtained 



for cooler stars, and C-enhancement or higher metallic ity are also enhanci ng t he F abundance 



or uppe r limits. Galactic chemical evolution models by 



RendaetaL 



(12004 and 



Kobavashi et al. 



(|2011a|) predict a similar F production at low metallicity, and are compared with observations. 
It is suggested that the v-process is necessary to explain the F production at the early stages of 
Galactic evolution, but the predicted [F/0] ratios d o not fit well the observed distribution, whereas 



the observed trend is closest to the SN II model of 



Renda et al. 



Kobayashi et ah (,2011a) and early phases of 



(|2004|) models with a low total energy in neutrinos, Ey = 3 x 10^^ e rg. Combined with 



2010|) with retrograde 



kinematics, HD 110281 is identified as a low-a star (iNissen & Schusteii i 
rotation. The consistent F nucleosynthesis of this star with the globular cluster halo stars support 
the scenario that this kind of objects may be accreted from nearby dwarf galaxies and enriched by 
early metal-deficient SN II. Future observations and modifications on models on relevant elements 
of the Galactic halo are important to achieve a better understanding of the above-mentioned results 
and remaining puzzles. 
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Fig. 1. — Overview of the HF-region spectra for all exposures of the programs stars. The position 
of the HF R9 line and the closeby C^^O are indicated by dash-dotted and dashed lines respectively. 
A K-hmd Sun spot spectrum in this area is also plotted for comparison. 
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Fig. 2. — - Na ab undances derived from CRIRES spectra are compared with those re-derived using 
Shetrone ( 1996)'s EWs measured from optical spectra (left), and the original results provided by 
Shetronel (|l996|) (right). Two objects with no EWs available from the literature are not included in 
the left plot. 
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Fig. 3. — Sample spectra with spectral synthesis of HF line at 23358 A of HD 110281 and 
HD 135148, the two program stars with clear HF detections. Black dots are observed data points. 
Solid lines refer to the synthetic spectra with the adopted A(F), with dashed lines and dotted lines 
referring to synthetic spectra with AA(F)=+ 0.1 dex. Note that for the CH star HD 135148, the 
C^^O line to the blue side of the HF line is clearly detected and also fitted. 



-26- 




Fig. 4. — Line bisectors of typical CO lines (solid lines) of different strength in the HF region. The 
line synthesis was done at constant C/0 ratio. Zero Doppler velocity corresponds to the laboratory 
wavelength. 
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Fig. 5. — Line bisectors of the HF R9 line (solid lines) of different strength. Zero Doppler velocity 
corresponds to the laboratory wavelength 
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Fig. 6. — Toomre diagram of the program stars. The long-dashed line corresponds to Vtotai = 
180 km s~^ The short-dashed line indicates zero rotation in the Milky Way. Two objects with HF 
detections, HD 1 10281 (with retrograde rotation) and HD 135 148 (the CH-star), and BD +01 2916 
(identified as a thick disk star) are marked in the plot. 
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Fig. 7. — Variances of the F abundance along 7^^, abundances of C, O and Fe. The filled circles 
with error bars refer to HD 110281 andHD 135148, and all upside-down filled triangles correspond 
to the program stars with 2-cr upper limits of F for the cases in which the HF line is not detected. 
1-cr errors are overplotted to the observed data. 
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Fig. 8. — The F yields in the form of [F/0] as a function of [0/H] and [0/Fe]. Models of 
Renda et al. (2004; thick lines) and Kobayashi et al. (2011a; thin lines) are shown, with line 
identifications described in the plot and te xt. Observations of vario us giants are als o inclu ded: 



open circles for field K and M giants from 



Cunha et al, 



(2003b, and 



Cunha & Smith 



amon ds, squares, crosses, and triangles fo r globular c 



M22 ( Alves-Britoetal. 



20121), NGC671 2 (lYong eta 



spectively; pluses for bulge giants from 



Cunha et al. 



uster giants from M4 (ISmith et alJ 



2008). and oj Cen (ICunha et al. 



(2005); di 



2005h . 



2003h . re 



(|2008h . Filled symbols are the same as in 



Figure |7l except that the CH-star HD 135148 is not included in the plot due to its peculiarly en- 
hanced F abundances as explained in the context. AH abundances of models a nd observations are 



calibrated to the same set of solar abundances from 



Anders & Grevessel ( 1989 ). 



Table 1: Stellar Parameters and Abundances of the Sample. 



Star 






[Fe/H] 


^km s ^ 


FWHM^ 


S/N 


A(C) 


A(0) 


A(F) 


A(Na) 


A(Mg) 


A(A1) 


A(Eu) 


BD +01 2916 


4200 


0.10 


-1.88 


2.10 


0.709 


357.3 


5.64 


7.11 


<2.19 


4.08 


6.21 


4.46 


-1.13 


HD 003008 


4100 


0.10 


-2.06 


1.90 


1.030 


591.8 


5.46 


7.09 


<1.67 


3.98 


6.33 


4.21 


-1.21 


HD 029574 


4200 


0.10 


-1.93 


1.75 


0.658 


226.8 


5.41 


7.18 


<2.34 


4.01 


6.30 


4.40 


-0.95 


HD 110281 


3950 


0.20 


-1.56 


2.43 


0.996 


462.5 


5.92 


7.76 


2.71 


4.33 


6.42 


4.70 


-0.46 


HD 118055 


4225 


0.70 


-1.78 


1.63 


0.687 


163.6 


5.77 


7.36 


<2.57 


4.03 


6.20 


4.58 


-0.86 


HD 135148 


4250 


0.70 


-1.96 


2.40 


0.899 


360.2 


6.96 


7.18 


2.59 


4.06 


6.02 


4.59 


b 


[S84] 2643 


4180 


0.40 


-2.13 


1.65 


0.735 


316.3 


5.37 


7.08 


<1.99 


4.09 


6.07 


4.50 


b 



"FWHM of the spectral lines of the sample stars (in mA), based on CO line measurements. 
''No Eu line or EW measurement is available for this object from 



Shetronel(ll996l) . 



Table 2: 3D Abundance Corrections^ on the HF R9 and Forbidden Oxygen Lines. 



Teff logg HFR9 [OI]6300A [O I] 6363 A 3D model'' 

3886 1.0 -0.02 -0.03 0.00 0.01 0.00 0.01 d3t36glOmm20n02 

4001 1.5 -0.03 -0.03 0.00 0.01 0.00 0.01 d3t40gl5mm20n01 

4478 2.5 -0.16 -0.12 -0.01 -0.01 -0.01 -0.01 d3t45g25mm20n02 

5024 2.5 -0.42 -0.52 0.00 -0.02 0.00 -0.02 d3t50g25mm20n01 



"Left columns give 3D - ID, right columns 3D - (3D), in dex. 
*AU models assume [M/H] = -2.0, and [or/Fe] = +0.4. 



